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Flavour Singlet Mesons in full QCD on the Lattice 
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We apply spectral methods to compute the OZI-rule suppressed loop-loop correlators in the pseudocalar meson 
flavour singlet channel. Using SESAM configurations (obtained with two degenerate sea quark flavours on 16 3 x 32 
lattices at /3 = 5.6, with standard Wilson action), we find for the first time clear evidence for mass plateau 
formation in the r\ channel of this theory. As a consequence, we observe a clear signal of a mass gap persistent 
under chiral extrapolation. This sets the stage for a more realistic two-channel approach, where partially quenched 
strange quarks would be included, in addition to the u, d sea quarks. 
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1. INTRODUCTION 

The large empirical mass gap between the rf 



and 7T masses, m,Q 



in. 



is commonly 



believed as being due to the Ua{^) anomaly of 
the flavour singlet axial vector current. In this 
picture, the large Ty'-meson mass is driven by 
quark-antiquark interactions mediated by topo- 
logical charge fluctuations in the QCD vac- 
uum. In the iV c = oo-limit of QCD this sce- 
nario is put forward through the famous Witten- 
Veneziano Q] formula (WVF) that models rriQ 
in terms of the topological susceptibility, from 
quenched QCD. The pitfalls of a lattice evalua- 
tion of the Witten-Veneziano formula have been 
elucidated in Ref. [Q. 

In principle an ah initio, direct verification of 
the Ua(1) origin of flavour singlet pseudoscalar 
masses beyond the large N c - approximation can be 
accomplished by simulating unqucnched QCD on 
the lattice, but conclusive results on the mass gap 
have been impeded to date ||], mainly due to se- 
vere signal-to-noise problems in the analysis of 
the flavour singlet correlator [||-[|. 

As shown in Fig. (fll), the correlator in question, 



CV(Ai) =< 0\O*(At)O(t = 0)|0 > 



(1) 



contains a loop-loop (fermionically 'discon- 
nected') contribution, D. This OZI-suppressed 
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Figure 1. The two contributions to the rf propa- 
gator in the full QCD vacuum: the connected part 
with running through valence quark lines and the 
'disconnected' piece where the valence quark lines 
annihilate instead of connecting source x and sink 
y at time slices t — and t = At. 



piece arises from the Wick contractions within 
each of the two light quark bilinears 



O = qj 5 q 



(2) 



and was evaluated in the literature by use of 
stochastic estimator techniques (SET), to cover 
its entire volume dependence. As a net result the 
77' propagator {At) is computed as the numer- 
ical difference between its connected and discon- 
nected pieces: 

CV(Ai) = C(At)-2D(At) A ^>°° exp(-rrvAt) , (3) 

and the ?7'-mass is to be extracted from its asymp- 
totic behaviour. Since the large mass m n i cor- 
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Figure 2. The net variation of low eigenmodes of 
Q over the sea quark mass range reached by the 
SESAM simulation, as expressed in the spread 
of eigenmode ratios, R = |A l (0.1575)/A i (0.1560)|, 
with the eigenvalues ordered according to increas- 
ing moduli. 



responds to a strong numerical cancellation be- 
tween the positive terms C and D, the calamity of 
previous LQCD attempts was to establish a win- 
dow in At with unambiguous single exponential 
decay, before signals disappear under the noise 
level. 

Barring brute force methods it appears that im- 
portant progress in the field could be achieved by 

1. avoiding SET and/or 

2. accomplish precocious ground state domi- 
nance in the flavour singlet channel. 

With this motivation, we report here on some 
recent methodological studies of ours where 
we rcaddrcss the unsettled problem of the direct 
T]' signal in LQCD and explore the possibility of 
spectral representations of quark propagators in 
superseding the SET treatment of D(At) @-|6). 
Such approach is of course in lign with the ex- 
pectation that chiral physics should be controlled 
by the low modes of the Dirac- Wilson operator. 
The point to find out is whether in the sea quark 
mass range of state-of-the-art full QCD simula- 
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Figure 3. The spectral cutoff dependence of the 
disconnected contribution, D l (At), for various 
time separations At = 1, • • • 16 (top to bottom), 
at the lightest SESAM sea quark mass. 



tions - like the SESAM simulation || - a rea- 
sonably small number of low eigenmodes actually 
suffices to saturate the fermion loop expansions 
with appropriate degree of accuracy. 

2. The spectral approach 

Since the Wilson-Dirac matrix, M, is non- 
normal at non-zero lattice spacings, it is prefer- 
able to expand in the eigenspectrum of its Her- 
mitian equivalent 

Qt = q .- 7sM . ( 4 ) 

Q has been shown to provide the optimal eigen- 
value basis in the sense of the singular value de- 
composition [fuj . 

For our practical benchmarking we use 200 
(195) SESAM lattices of size 16 3 x 32 with k = 
0.1575 (0.1560) |,||. A detailed account of our 
analysis can be found in Ref. || . 

The eigenmodes of Q can be determined by 
means of Lanczos methods. In Fig. (||) the 300 
lowest eigenvalues are plotted to show their re- 
sponse under the SESAM variation in sea quark 
mass. The characteristic range of this response in 
eigenvalue space suggests that a truncated eigen- 
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Figure 4. The disconnected contribution to the 
^'-correlator from its spectral representation with 
cutoff I — 300, at the smallest sea quark mass of 
the SESAM simulation. 



mode approach (TEA), based on some few hun- 
dred low modes, might indeed encompass the in- 
frared physics of interest. 

In terms of the eigenstates of Q, 

Qm = \i\ipi), (5) 

the quark propagator can readily be computed, 
from any source z to every sink z'\ 

M -1 (z,z') = 

^ 75 |^(z))(^(zQ| 

From this expression it is straightforward to com- 
pute the loop-loop correlator of Fig. ([!]) in spec- 
tral representation, T; it should be dominated by 
a small number of modes at large enough values of 
At. This is indeed reflected in Fig. (f|) where we 
display the Z-dependence of the disconnected con- 
tribution to the ry'-propagator, D l , the family of 
curves representing the set of different time slices, 
At = 1, ■ ■ • , 16 (from top to bottom). While we 
find a few tens of modes to bear out the long 
range (in At) features of D, a few hundred modes 
appear to saturate its short range structure un- 
der the SESAM conditions. This finding is also 



Figure 5. The analogue to Fig. (3) , but for the 
connected contribution, C(At), to Eq. (3). 



visualized in Fig. (^) where we plotted D 300 (At) 
together with the previous SET result Q and find 
them to agree nicely within their errors. 

For the connected piece of the r/-correlator, C, 
the viability of TEA is less evident as can be seen 
from Fig. (||). Even with a cutoff I = 300, sat- 
uration at large At is not fulfilled in detail, not 
to speak of the region of small time separations, 
where excited state contributions are not covered 
at all by the 300 lowest modes. This is of course 
not an obstacle for our present purposes, since 
this part of the correlator can easily be treated by 
standard methods. In fact, at this point, we go 
one step beyond and exploit the high accuracy of 
C as achieved from iterative solvers: we replace C 
by its ground state contribution, (as obtained 
from a single exponential fit at large time separa- 
tions, see Fig. (^)), C — > G%. This then leads us 
to the final form of the projected ?/-correlator 

C n ,(At) = C3(At)-2T{At) . (7) 

in terms of the two loop spectral approximation, 
T: 

T(At)=tr\ t=0 Q- 1 tr| At Q-\ (8) 

with Nf — 2 being the number of active sea quark 
flavours and Q~ x taken from Eq. (||). 
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Figure 6. The one loop correlator computed by 
standard linear solver techniques (labeled 'stan- 
dard method') and TEA with spectral cutoff I = 
300, at the lightest SESAM sea quark mass. Also 
shown is the ground state projection (marked 
'Fit'). 



Figure 7. Effective rj'-m&ss in lattice units at the 
lightest (k = 0.1575) and heaviest (k = 0.1560) 
sea quark masses of the SESAM simulation, both 
with (SS) and without (LL) source and sink 
smearing 



This enables us to carry out a standard effec- 
tive local ?/-mass analysis, which is much more 
sensitive than looking at correlators or ratios 
thereof Q . As a result we find an effective mass 
plot (see Fig. ([?])) with striking mass plateau 
formation which extends over the entire interval 
1 < At < 9, for the smallest sea quark mass. 
This precocious ground state dominance in the 
T)' channel provides us with hitherto unobserved 
accuracy for the 7y'-mass estimate, over the entire 
range of sea quark masses of the SESAM QCD 
simulation! 

Note from Fig. (Q) that the flavour-nonsinglet 
mass (marked V) is clearly resolved from the rj- 
mass plateau, with increasing mass gap as the 
quark mass is lowered, at our value of lattice spac- 
ing, a- 1 {m p ) = 2.30 GeV. 

Since TEA works better with decreasing sea 
quark masses and Lanczos methods do not dete- 
riorate with growing condition numbers, TEA is 
expected to be by far superior to SET in the deep 
chiral regime of tomorrow's simulations where lin- 
ear solvers will tend to become unstable. 



3. Discussion and outlook 

Our analysis of the rj' signal in two-flavour 
QCD has led us to the reassuring result that spec- 
tral methods in conjunction with ground state 
projection on the connected contribution of the 
//-correlator can resolve the rj'-w mass gap be- 
fore the chiral and continuum extrapolations. In 
Fig. (H) we show that the lattice data are suf- 
ficient to sustain a chiral extrapolation as well. 
The figure demonstrates moreover that there is 
full consistency of TEA with SET data when be- 
ing analysed in the same manner. Before con- 
tinuum extrapolation, a linear fit of rn^, (m sea 
works with reasonable % 2 and results in M n > = 
307(15) MeV. This value is far away from the ac- 
tual experimental number, m v ' = 957 MeV; this 
does not come as a surprise, however, since the 
SESAM QCD vacuum configurations are miss- 
ing the strange sea quark contributions. There- 
fore our pseudoscalar flavour singlet meson resem- 
bles more an (octet) rj than the rj meson proper. 
Apart from the continuum extrapolation t.b.d., it 
is therefore necessary to introduce strange quarks 
into the calculation, in order to make contact with 
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Figure 8. Demonstration of the quality of data: 
chiral extrapolation of rn^,, plotted vs. sea quark 
mass msea, all in lattice units. m pe is a naive 
estimate from WVF (replacing m§ — * 2/3mQ to 
accomodate to Nf = 2). 



real experiment. 

Since Nf = 3 simulations are not feasible at 
this time the next step would be to consider both 
r\ and rf in a partially quenched setting with two 
degenerate light quarks n in the sea, and a strange 
valence quarks added. This opens the scenario of 
a coupled two channel approach of pseudoscalar 
singlet mesons in the n and s sectors, in terms of 
the quark flavour basis 

|r?> = cos$ l^nn > +sin$ |?7 SS > (9) 
\rf > = -sin$ |77nn > +cos$ |77 SS > 

with Eq. (||) being replaced by a 2 x 2 correlator 

C(i\ = ( Gnn{ $ ~ 2Z3 ™«W -V2D ns (t) 
U ^ -V2D sn (t) C„(t) - D ss (t) 

So far we have dealt only with the left hand up- 
per corner of this matrix. The evaluation of the 
OZI-rule violating contributions in its remaining 
entries are expected to proceed along the lines 
described above, without any additional compli- 
cation. The mixing problem then amounts to the 
solution of the eigenvalue problem |l3| 

C{t)C- 1 (t )\i>= X l {At)\i> withi>£ (10) 



with t — to = At. The two eigenmodes coincide 
with the physical rj' and rj states. Their masses 
are to be extracted through the asymptotic rela- 
tions 

A; (At) A -^°° exp(-m 2 At) . (11) 

Previously, the 77-77' mixing could only be ad- 
dressed by perturbative modeling j| . One would 
expect that the precocious mass plateau forma- 
tion established in the Tin-sector will persist in the 
presence of channel mixing from this eigenvalue 
problem. Work along this line is in progress jl4| . 
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